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A  generic,  transient  fuel  cell  ohmic  loss  mathematical  model  was  developed  for  the  case  of  contaminants 
that  ion  exchange  with  ionomer  protons.  The  model  was  derived  using  step  changes  in  contaminant 
concentration,  constant  operating  conditions  and  foreign  cation  transport  via  liquid  water  droplets.  In 
addition,  the  effect  of  ionomer  cations  redistribution  within  the  ionomer  on  thermodynamic,  kinetic  and 
mass  transport  losses  and  migration  were  neglected.  Thus,  a  simpler,  ideal,  ohmic  loss  case  is  defined 
and  is  applicable  to  uncharged  contaminant  species  and  gas  phase  contaminants.  The  closed  form  solu¬ 
tions  were  validated  using  contamination  data  from  a  membrane  exposed  to  NH3.  The  model  needs  to 
be  validated  against  contamination  and  recovery  data  sets  including  an  NH4+  contaminated  membrane 
exposed  to  a  water  stream.  A  method  is  proposed  to  determine  model  parameters  and  relies  on  the 
prior  knowledge  of  the  initial  ionomer  resistivity.  The  model  expands  the  number  of  previously  derived 
cases.  Most  models  in  this  inventory,  derived  with  the  assumption  that  the  reactant  is  absent,  lead  to 
different  dimensionless  current  vs.  time  behaviors  similar  to  a  fingerprint.  These  model  characteristics 
facilitate  contaminant  mechanism  identification.  Separation  between  membrane  and  catalyst  (electroin¬ 
active  contaminant)  contamination  is  conceivably  possible  using  additional  indicative  cell  resistance 
measurements.  Contamination  is  predicted  to  be  significantly  more  severe  under  low  relative  humidity 
conditions. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Proton  exchange  membrane  fuel  cells  (PEMFCs),  an  alternative 
power  source  currently  being  developed  [1],  is  sensitive  to  con¬ 
taminants  present  in  reactant  streams  and  released  from  system 
components  [2].  Generic  models  were  developed  to  address  the 
large  number  of  existing  and  unidentified  contaminants  [3-6].  Such 
models  are  useful  to  classify  contaminants  (mechanism  identifica¬ 
tion),  set  tolerance  limits  and  design  experimental  procedures  for 
validation  and  extraction  of  model  parameters.  For  instance,  fuel 
cell  testing  with  low  contaminant  concentrations  requires  signif¬ 
icant  time  and  resources  because  effects  are  small  and  not  easily 
measured  with  certainty.  A  simple  model  for  membrane  contam¬ 
ination  has  not  yet  been  derived  and  is  needed  to  complete  the 
existing  library. 

Several  models  have  already  been  proposed  showing  that  all 
types  of  performance  losses  take  place  within  a  contaminated 
PEMFC  including  thermodynamic,  kinetic,  ohmic  and  mass  trans¬ 
port  losses  [7-13].  Experimental  data  support  this  statement  [14]. 
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However,  in  most  model  cases,  the  physics  associated  with  the  ion 
exchange  process  were  not  taken  into  account  [7-12].  Langmuir 
isotherms  were  recently  used  [13]  but  such  an  approach  was  dis¬ 
continued  to  characterize  ion  exchangers  [15]  because  it  does  not 
reproduce  their  behavior. 

A  simple  transient  PEMFC  contamination  model  is  derived  using 
a  similar  approach  as  prior  efforts  [3-6]  with  a  focus  on  the  mem¬ 
brane  ion  exchange  with  foreign  cations.  The  approach  considers 
first  order  effects,  to  maintain  model  generality  and  applicability 
to  many  different  cations,  and  the  derivation  of  an  analytical  solu¬ 
tion,  to  simplify  its  use  for  predictions,  validation  and  extraction  of 
model  parameters. 

2.  Model  assumptions 

Fig.  la  illustrates  the  path  taken  by  a  foreign  cation  X+n  from 
the  flow  field  channel  to  the  ionomer  exchange  site.  It  is  assumed 
that  liquid  water  droplets  maintain  a  steady  flow  and  concentra¬ 
tion  of  foreign  cations  originating  from  the  environment  or  system 
components.  Liquid  water  is  not  necessarily  required  to  introduce 
foreign  cations  into  the  fuel  cell.  Particles  containing  soluble  salts 
and  entrained  by  the  reactant  flow  provide  a  supply  of  foreign 
cations  by  dissolution  in  fuel  cell  areas  containing  water  droplets. 
The  presence  of  product  water  in  the  gas  diffusion  electrode  creates 
a  path  for  the  foreign  cations  to  reach  the  ionomer  by  diffusion.  The 
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Nomenclature 

a  empirical  Eqs.  (16)  and  (17)  parameter  (respectively 

cm-1  and  dimensionless) 

A  membrane  surface  area  (cm2) 

b  empirical  Eqs.  (16)  and  (17)  parameter  (respectively 

cm-1  and  dimensionless)  orTafel  slope  (V) 
cg  foreign  contaminant  concentration  in  the  gas  phase 

(mol  cm-3) 

Cj  equivalent  ionic  concentration  in  the  ionomer  phase 

(mol  cm-3) 

C[  h+  proton  concentration  in  the  ionomer  phase 
(mol  cm-3) 

ci  x+n  foreign  cation  concentration  in  the  ionomer  phase 
(mol  cm-3) 

q  equivalent  ionic  concentration  in  the  liquid  phase 

(mol  cm-3) 

q  H+  proton  concentration  in  the  liquid  phase  (mol  cm-3 ) 

C/  x+n  foreign  cation  concentration  in  the  liquid  phase 

(mol  cm-3) 

Cr  reactant  concentration  (mol  cm-3 ) 

d  distance  (cm) 

dg  gas  phase  boundary  layer  thickness  (cm) 

di  ionomer  phase  thickness  (cm) 

d/  liquid  phase  thickness  (cm) 

Dg  diffusion  coefficient  in  the  gas  phase  (cm2  s-1 ) 

Df  diffusion  coefficient  in  the  ionomer  phase  (cm2  s-1 ) 

D/  diffusion  coefficient  in  the  liquid  phase  (cm2  s-1 ) 

e  polarization  curve  parameter  characterizing  mass 

transport  effects  (V) 

E  cell  potential  (V) 

Ea  anode  potential  (V) 

Ec  cathode  potential  (V) 

Eq  polarization  curve  constant  cell  potential  parameter 

(V) 

Estm  experimental  estimates  for  specific  model  parame¬ 

ters  (m  is  an  integer  with  a  value  of  1-4) 

/  polarization  curve  parameter  characterizing  mass 

transport  effects  (cm2  A-1 ) 
i  current  density  (A  cm-2 ) 

ix  current  density  associated  with  a  contaminant  X 

(A  cm-2) 

ix+n  current  density  associated  with  a  contaminant  X+n 
(A  cm-2) 

n  foreign  cation  charge 

N  proton  flux  in  the  liquid  phase  (mol  s-1  cm-2) 

Pe  Peclet  number 

R  areal  membrane  resistance  ( £2  cm2 ) 

t  time  (s) 

t0. 5  time  required  for  the  dimensionless  current  to 

decrease  to  a  value  equal  to  half  the  difference 
between  initial  and  steady  state  values  (s) 
u  liquid  phase  velocity  (cm  s-1 ) 

V  ionomer  volume  (cm3 ) 

x  equivalent  ion  fraction  in  the  liquid  phase 

xH+  equivalent  proton  fraction  in  the  liquid  phase 
xx+n  equivalent  foreign  cation  fraction  in  the  liquid  phase 

y  equivalent  ion  fraction  in  the  ionomer  phase 

yH+  equivalent  proton  fraction  in  the  ionomer  phase 
yx+n  equivalent  foreign  cation  fraction  in  the  ionomer 

phase 

a^tn  separation  factor 

a  ionomer  conductivity  ( cm-1 ) 


a 


0  d,  df+dj 


Flow  Gas 

field  diffusion 

channel  electrode  Membrane 


Fig.  1.  Contaminant  X+n  transport  pathway  towards  the  ionomer  (a).  Liquid  phase 
transport  in  the  flow  field  and  the  gas  diffusion  electrode  is  followed  by  pene¬ 
tration  into  the  ionomer  at  its  surface.  Subsequently,  X+n  diffuses  towards  an  ion 
exchange  site  where  it  associates.  Ion  exchange  stoichiometry  is  not  taken  into 
account  in  this  schematic  illustration.  Equivalent  ion  fractions  in  both  liquid  and 
ionomer  phases  under  liquid  phase  mass  transport  control  after  contamination  ini¬ 
tiation  and  before  a  steady  state  is  reached  (b).  For  the  case  illustrated,  the  foreign 
cation  X+n  is  preferred  by  the  ionomer  <  1  and  yx+n  >  xx+n). 


existence  of  such  paths  at  preferential  locations  on  the  gas  diffu¬ 
sion  electrode  was  observed  with  in  situ  transparent  cell  designs 
and  ex  situ  [16,17].  The  path  length  is  assumed  to  be  of  the  same 
order  as  the  electrode  thickness  (Fig.  la).  At  the  water/ionomer 
interface,  the  foreign  cations  easily  enter  the  ionomer.  This  process 
is  not  rate  determining  for  theoretical  reasons  and  is  supported  by 
experimental  data  [15].  Subsequently,  foreign  cations  move  within 
the  ionomer  bulk  by  diffusion  before  associating  themselves  with 
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Fig.  2.  Schematic  ion  exchange  isotherms  for  the  case  of  a  foreign  cation  X+n 
and  H+  (Nation  in  contact  with  a  dilute  solution  of  XYn).  Accessible  xx+n  range 
during  contamination  and  recovery  are  also  indicated.  An  isotherm  of  the  form 
yx+n  =  1  /[  1  +aj+n(l  -  xx+n)/xx+n ]  was  used  assuming  equal  cation  charges  and  a 
constant  selectivity  coefficient  thus  reducing  the  expression  to  the  separation  fac¬ 
tor  definition  [15,32].  For  illustrative  purposes,  a^+n  =  1.2  for  Li+  whereas  for  other 
cations,  a^_n  =  0.4. 


exchange  sites.  This  last  step  is  not  rate  determining  unless  the 
exchange  sites  chelate  foreign  cations  to  form  complexes  [15].  The 
protons  leaving  the  ionomer  travel  through  the  entire  electrode 
thickness  to  compensate  for  the  foreign  cations  flux  towards  the 
ionomer  and  maintain  the  charge  balance  (the  anion  concentration 
is  constant,  Fig.  lb). 

Anions  were  not  included  in  Fig.  la  because  they  are  prevented 
from  entering  the  ionomer  owing  to  Donnan  exclusion  [15,18].  The 
presence  of  fixed  negative  charges  (sulfonate  sites)  within  the  poly¬ 
mer  ensures  the  absence  of  foreign  anions  for  solutions  with  a 
concentration  up  to  0.1  N  or  more.  Water  production  within  fuel 
cells  ensures  that  in  most  cases  this  limit  is  not  exceeded  unless 
a  catastrophic  event  takes  place  (for  example,  flooding  of  the  air 
intake  by  sea  water).  Therefore,  it  is  assumed  that  the  concentra¬ 
tion  of  the  foreign  cation  is  relatively  low  (anions  do  not  need  to 
be  considered).  By  contrast,  Donnan  exclusion  is  not  always  con¬ 
sidered  to  explain  the  anion  effect  [19].  This  statement  does  not 
imply  that  anions  do  not  have  any  effect.  For  instance,  the  ionomer 
proton  spills  over  the  catalyst  surface  and  diffuses  on  its  surface 
thus  enlarging  the  active  area  outside  the  ionomer  [20-22].  The 
presence  of  a  complexing  anion  dissolving  the  Pt  would  reduce  the 
active  area  and  increase  the  kinetic  losses.  Such  kinetic  losses  are 
not  considered  here. 

Foreign  cations  such  as  NFf4+,  Na+,  K+,  Rb+,  Cs+,  Ag+,  Ca2+,  Co2+, 
Cu2+,  Ni2+,  Zn2+,  Fe3+  and  Ce3+  are  generally  preferred  by  the  Nation 
ionomer  over  FT  [12,23-31].  Only  Li+  appears  to  show  a  differ¬ 
ent  behavior  (FT  is  slightly  preferred  [12,27])  but  this  observation 
does  not  affect  membrane  contamination  deductions.  Representa¬ 
tive  ion  exchange  isotherms  for  these  specific  cases  are  depicted  in 
Fig.  2  to  discuss  the  extent  of  the  foreign  cation  exchange  with  the 
ionomer  proton.  Other  ions  also  need  to  be  considered  in  the  prod¬ 
uct  water  including  the  foreign  anion  Y-  and  water  dissociation 
products: 

2H20  H30+  +  0H-  (1) 


It  is  assumed  that  the  foreign  anion  is  not  hydroxyl  and,  X(OH)n 
and  FIY  are  highly  soluble.  The  equivalent  foreign  cation  and  proton 
ionic  fractions  are  [15]: 

ncix+n  ncux+n 

Xx+n  =  ^ - =  — 7 —  W 

nci,x+n  +  c/,h+  ci 


XH+  - 


Cf,H+ 

ncl,X+n  +  Cj  H+ 


1  —  Xx+n 


(3) 


Fig.  3.  Validation  of  the  selected  ion  exchange  isotherm  ( yx+n  =  1/[1  +aj^n(l  - 
Xx+n)/xx+n])  with  different  cations  in  Nation  417  data  [31].  Separation  factors  and 
correlation  coefficients  are  given  in  Table  1. 


During  the  contamination  process,  if  nqx+n  >>  clH+  then 
Xx+n  &  1  and y x+n  «  1  whereas  if  nqx+n  <<  clH+  thenxx+n  «  Oand 
yx+n  &  0  (Fig.  2).  For  other  intermediate  values  of  xx+n,  generally 
yx+n  >  xx+n  (Fig.  2).  In  other  words,  it  is  assumed  that  the  contami¬ 
nant  concentration  may  span  the  entire  range  (Fig.  2).  Therefore,  to 
minimize  the  membrane  contamination  extent,  the  foreign  cation 
concentration  needs  to  be  significantly  lower  than  the  proton  con¬ 
centration  in  the  product  water  (10_7M  at  25  °C).  The  situation 
during  the  recovery  process  is  different  because  nqx+n  <<  q  H+, 
Xx+n  ~  0  and  yx+n  «  0  (Fig.  2).  In  other  words,  few  foreign  cations 
are  present  near  the  ionomer.  Thus,  it  is  expected  that  a  contami¬ 
nated  cell  will  recover  most  of  its  performance  loss  without  adding 
any  mitigation  strategy  if  given  sufficient  time  and  contact  with 
liquid  water.  The  present  discussion  is  extendable  to  other  cases 
including  other  foreign  anions  (different  species  and  ionic  charges), 
existence  of  a  common  anion  (OF!-),  sparsely  soluble  components 
and  the  presence  of  C02  in  air  (pFI  change  affecting  Eqs.  (2)  and 
(3)): 

C02  +  2H20  HC03-  +  H30+  (4) 


The  ion  exchange  isotherm  is  expressed  by  the  separation  factor 
[15]: 

H+  _  yH+*x+”  _  (i  -yx+n)xx+n 
x+n  yx+"*H+  yx+n(i-*x+n) 

Eq.  (5)  leads  to: 

yx+n  =  - — ,  h+  m  1 - w — ; 

Eq.  (6)  should  not  be  confused  with  the  selectivity  coefficient 
which  is  similar  in  form  to  a  chemical  equilibrium  constant  and 
includes  stoichiometric  coefficients  arising  in  the  ion  exchange 
equilibrium  due  to  unequal  cation  charges  [15].  Eq.  (6)  corresponds 
to  the  selectivity  coefficient  assuming  equal  cation  charges  [15,32]. 
The  selectivity  coefficient  only  mimics  the  chemical  equilibrium 
constant  because  ion  exchange  is  not  a  chemical  reaction  (ions  do 
not  react  but  change  phases).  As  a  result,  the  selectivity  coefficient 
is  not  constant  [15]  as  opposed  to  a  chemical  equilibrium  constant. 
Fig.  3  data  indicate  that  Eq.  (6)  ion  exchange  isotherm  accurately 
represents  the  behavior  of  an  ionomer  in  the  presence  of  different 
cations  even  if  in  several  cases  cations  do  not  have  equal  charges. 
Adoption  of  Eq.  (6)  contributes  to  the  development  of  a  generic 
contamination  model. 

Foreign  cations  move  in  the  product  water  by  two  different 
migration  effects.  The  product  water  moving  in  gas  diffusion  elec¬ 
trode  pores  (Fig.  la)  is  subjected  to  a  potential  gradient  extending 
from  the  catalyst  layer  to  the  flow  field  channel.  This  potential  gra- 
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/  /  A  cm  2 

Fig.  4.  Schematic  polarization  curves  during  the  contamination  process  at  constant 
cell  potential.  The  operating  points  are  located  in  the  ohmic  regime.  The  initial 
operating  point  gradually  moves  from  point  A  to  point  B  as  the  ionomer  is  con¬ 
taminated  with  the  foreign  cation  and  the  resistance  increases.  The  polarization 
curves  were  generated  using  E0-b  log (z)  -Ri-e  exp (fi)  with  E0  =  0.95  V,  b  =  0.07  V, 
e  =  0.1  V,/=  1  cm2  A-1  ( R  values  are  given  in  figure). 


dient  is  small  and  typically  less  than  1 0-20  mV  at  ~1  A  cm-2  [33,34] 
and  thus  is  not  considered  an  important  model  element  because  the 
water  ionic  resistivity  is  high  and  the  resulting  current  is  small.  The 
contribution  arising  from  the  difference  in  cation  mobility  is  not 
negligible  (H+,  X+n)  and  is  assessable  using  liquid  junction  poten¬ 
tial  theory  [35].  Cation  migration  in  dilute  solutions  such  as  ground 
water  is  negligible  because  the  contaminant  ionic  strength  is  much 
smaller  than  other  ions  [36]  acting  as  a  supporting  electrolyte  [35]. 
However,  this  assumption  cannot  be  made  here  otherwise  it  would 
minimize  the  extent  of  contamination  (Eq.  (6))  by  reducing  the 
foreign  cation  equivalent  ionic  fraction  (Eq.  (2)).  Ion  migration  is 
assumed  to  be  negligible  here  for  several  reasons  including  consid¬ 
eration  of  relevant  cases  such  as  uncharged  contaminant  uptake 
by  the  ionomer  affecting  its  conductivity,  gas  phase  contaminants 
(Section  4.1)  and  accuracy  determination  of  a  simpler  model.  As 
already  mentioned,  the  main  objective  is  the  development  of  an 
explicit  analytical  solution  to  the  contamination  model.  Such  an 
objective  is  not  considered  possible  with  the  inclusion  of  migra¬ 
tion  effects.  For  instance,  an  analytical  solution  to  the  ion  exchange 
problem  is  only  possible  for  two  cases  [15]  requiring  restrictive 
assumptions  such  as  absence  of  ionomer  selectivity  (a^+n  =  1 )  and 
equal  valence  ions  or  equal  counter  ions  mobility.  Convection  is  also 
assumed  to  be  negligible  because  the  Peclet  number  Pe  =  udi/Di  [37] 
is  less  than  1  (u  <  0.0005  cm  s-1  at  0.1  A  cm-2  [38],  d/  =  0.02  cm,  and 
D/  =  10-5  cm2  s-1). 

The  current  density  is  assumed  to  be  relatively  low  to  minimize 
cation  distribution  gradients  within  the  ionomer  causing  thermo¬ 
dynamic,  kinetic,  and  mass  transport  losses  [7-13].  These  gradients 
are  also  expected  to  affect  ohmic  losses  as  well  as  the  ion  exchange 
process  because  at  the  ionomer  interface  the  foreign  cation  concen¬ 
tration  is  larger  than  the  average  value  (water/ionomer  boundary 
condition).  The  simplification  introduced  (relatively  low  current 
density)  ensures  that  all  these  effects  are  minimized  thus  defining 
an  ideal  ohmic  loss  case  useful  to  study  the  more  general  case  and 
applicable  to  a  few  specific  cases  (uncharged  contaminant  uptake, 
gas  phase  contaminants).  It  is  also  assumed  that  the  cell  potential 
is  controlled  within  the  ohmic  regime  (Fig.  4).  The  cell  polariza¬ 
tion  curve  undergoes  modifications  as  the  ionomer  is  contaminated 
by  the  foreign  cation.  The  behavior  illustrated  in  Fig.  4  reproduces 
experimental  data  for  different  ionomer  contamination  levels  [14]. 
As  contamination  progresses,  the  cell  operating  point  moves  from 
initial  point  A  to  point  B  while  the  trajectory  remains  in  the  ohmic 
regime. 

Analysis  of  all  transport  processes  indicate  that  transport  only 
occurs  by  diffusion  in  both  phases.  Additionally,  only  liquid  phase 


or  ionomer  phase  diffusive  transport  is  rate  determining.  A  quan¬ 
titative  criteria  was  derived  to  determine  the  rate  controlling  step 
based  on  the  following  dimensionless  group  [15]: 


CiDidi 


<<  1,  ionomer  phase  diffusion  control  (7) 


(5  +  2oix+n)  >>  1,  liquid  phase  diffusion  control  (8) 


C/D/dj 


If  the  diffusion  coefficients  are  significantly  different,  a  harmonic 
mean  can  be  used  [15].  Here,  q  »  q  with  q  =  nci  x+n  +  ciH+  (dilute 
contaminant  of  >10-7  M  level  to  ensure  ion  exchange  process  takes 
place  and  low  ionomer  equivalent  weight  leading  to  an  exchange 
site  concentration  of  M  level),  Dl  (well  humidified  membrane 
with  expanded  pores),  d/  <$c  d;  (25-50  [xm  membrane  thickness 
in  comparison  to  ~200p,m  gas  diffusion  electrode  thickness  as 
depicted  in  Fig.  la)  and  a^+n  <  1  (for  most  foreign  cations  with 
the  exception  of  o£!+~l).  Therefore,  contamination  by  a  foreign 
cation  is  liquid  phase  controlled  (Eq.  (8)  applies,  Fig.  lb)  because 
qDjdi/qDidi  1  by  several  orders  of  magnitude  and  5  +  2o^+n  >  1 . 

Operating  conditions  are  assumed  to  be  either  kept  constant 
(stoichiometry,  relative  humidity,  cell  potential)  or  controlled 
within  a  narrow  range  (temperature,  pressure  [39])  to  maintain 
a  constant  water  balance  [40]  and  transport  parameter  values  (dif¬ 
fusion  coefficient).  It  is  also  assumed  that  the  chemical  membrane 
degradation  affecting  its  ion  exchange  capacity  (loss  of  exchange 
sites)  takes  place  with  a  longer  time  scale  than  cation  contamina¬ 
tion  [41,42]  which  is  favored  by  the  use  of  wetter  reactant  streams 
[43].  This  statement  is  important  considering  that  foreign  cations 
such  as  Cu2+,  Fe2+  and  Fe3+  accelerate  thermal  [44]  and  chemical 
[45-49]  ionomer  degradation. 


3.  Model  equations 


An  ion  exchange  model  for  diffusion  control  in  the  liquid  phase 
(absence  of  migration  and  convection)  was  derived  [  1 5  ].  This  model 
is  reformulated  here  for  the  case  of  an  ionomer  sheet.  Fields  first  law 
for  the  proton  flux  is: 


N  =  D, 


c,  H+(d,,t)-c/jH+( o,  t) 

d, 


(9) 


The  membrane  proton  mass  balance  is: 


dcUH+(t)V 

dt 


=  -AN 


(10) 


Reformulation  of  Eq.  (6)  in  terms  of  concentrations  leads  to  the 
following  equilibrium  condition  at  the  liquid  phase/ionomer  phase 
interface: 


ci,  h+(0 


Q«X+nc|, H+(dh  t) 

Cl  +Cl,H+(dh  t)(0#+n  -1) 


(11) 


It  is  assumed  that  the  ionomer  is  initially  in  a  fully  protonated 
form  and  that  the  foreign  cation  is  much  more  concentrated  than 
the  proton  concentration  in  the  flow  field  channel  water: 


c,H+(d,,0)  =  c, 
ci,H+(9>  t)  =  0 

Combination  of  Eqs.  (9)  to  (1 1 )  and  (13)  leads  to: 

dc,  H+(d),  t)  / _ l _ 

dt  \cl  H+(dh  t)(c,  +  (<*”+„  -  l)c,  H+(d,,  t)) 

_ («x+n  -i)q,H+(rf),t) _ j  _  — D| 

c;,H+(d|,  t)(c,  +  (afn  -  1  )c,  H+(d,,  t))2  /  Cjdjd/a^ „ 


(12) 

(13) 


(14) 
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The  solution  to  Eq.  (14)  with  the  boundary  condition  given  by 
Eq.  (12)  and  re-use  of  Eq.  (11)  to  replace  q  H+(d/ ,  t)  by  q  H+(t)  is: 


ln(yH+(t))  + 


1 

q;H+ 

«X+" 


(i  -yH+(t))  = 


-qDi 

Cidid,afn 


(15) 


The  ionomer  cation  composition  is  linked  to  its  conductivity  by 
a  linear  (NH4\  Na+,  Ca2+  [23,25,28,30,31 ,50])  or  logarithmic  (NH4+, 
Na+,  Ca2+,  Fe3+  [30,31])  empirical  relationship: 


cr  =  a  +  byH+  (16) 

In  o  =  a  +  byH+  (17) 


Eq.  (16)  only  applies  to  the  absence  of  interactions  between  the 
proton  and  the  foreign  cation.  Other  relationships  may  be  used 
because  Eqs.  (16)  and  (17)  do  not  appear  to  fit  Cu2+,  Ni2+  and  Fe3+ 
data  [26].  For  demonstration  purposes,  combination  of  Eq.  ( 1 6)  with 
Eq.  (15)  leads  to: 


cr{t)-a\ 

b  J 


~ciDi  t 

Cjdjd,afn 


(18) 


The  cell  voltage  (controlled)  and  areal  membrane  resistance  are, 
respectively: 


E  =  EC-Ea-Ri  (19) 

R= ^  (20) 

<7 

Combination  of  Eqs.  (19)  and  (20)  leads  to  (changes  in  ther¬ 
modynamic,  kinetic  and  mass  transport  losses  are  assumed  to  be 
negligible): 

i(t )  _EC-Ea-E  R( 0)  _  R(0)  _  dfor(t)  _  a(t) 

i(  0)  R(t )  Ec-Ea-E  R(t)  G{0)di  a{  0)  1  J 

A  similar  reduction  is  not  possible  for  the  galvanostatic  case. 
Furthermore,  Eq.  (21 )  is  consistent  with  prior  generic  models  [3-6]. 
Combination  of  Eqs.  (16),  (18)  and  (21)  leads  to: 


ln^(i(t)/i(0))(a  +  b)-aj 


+ 


(i(t)/i(0))(a  +  b)-a\ 

b  J 


-QPj 

Cidid,ax+n 


-t  (22) 


For  0  and  t-y  oo,  Eq.  (22)  reduces  to,  respectively: 

Kt)  _  1 _ b  c,D,  t 

i(0)  a  +  b  Cjd/d/ 

ijt)  _  a 

i(0)  a  +  b 


(23) 

(24) 


For  the  recovery  situation,  it  is  assumed  that  the  ionomer  is  ini¬ 
tially  in  a  fully  foreign  cation  form  and  that  the  proton  is  much  more 
concentrated  than  the  foreign  cation  concentration  in  the  flow  field 
channel  water: 


c/,H+(d/>  0)  =  0  (25) 

0,h+(6>  0  =  Q  (26) 

Combination  of  Eqs.  (9)  to  (11)  and  (26)  leads  to: 


dcf  H+(d/,  t)  / _ 1 _ 

dt  v(q,H+(di  -  0  -  q)(q  +  («£+„  -  i)q,H+(d/,  9) 

_ («^+„  -l)cliH+(d|,t) _ \  _  -D; 

(c,  H  Ad,,  t )  -  c,Xc,  +  (af„  -  1  )e,  H+(dh  t)f  J  c,d,d,  af+n 

(27) 


Fig.  5.  Effect  of  contamination  by  a  foreign  cation  and  subsequent  recovery 
on  the  cell  normalized  current  density  (Eqs.  (22)  and  (29)).  a  =  0.01  £2-1  cm-1, 
b  =  0.09£2_1  cm-1,aj+n  =0.1,  Cidjdi/c/D/  =  lOh. 


The  solution  to  Eq.  (27)  with  the  boundary  condition  given  by 
Eq.  (25)  and  re-use  of  Eq.  (11)  to  replace  q  H+(d/,  t)  by  cf  H+(t)  is: 

In (1  -yH+(t))+  (l  -  t^x+n)  yH+(f)  =  f  (28) 

Combination  of  Eqs.  (16)  and  (21)  with  Eq.  (28)  leads  to: 
(i(t)/i(0))(a  +  b)-cA 

J 


In  1  - 


,  f ,  „,H+  ,  ( (i(t)/i(0 ))(a  +  b)-a\  _  -qD,a^n 

+  [1  ~ax+n)  (  b  J  ~  Cjdjd/ 


(29) 


For  t^O  and  oo,  Eq.  (29)  reduces  to,  respectively: 

121  -  a  ,  b  ci°i  t 

i(0)  a  +  b  a  +  b  qdjd/ 


(30) 

(31) 


Fig.  5  illustrates  Eqs.  (22)  and  (29)  for  arbitrary  parameter  values. 
It  is  observed  that  recovery  is  significantly  slower  than  contami¬ 
nation,  which  is  consistent  with  a  preferred  foreign  cation  in  the 
ionomer  (Figs.  2  and  3),  but  proceeds  to  completion  with  the  pro¬ 
ton  form  re-established  because  xx+n  ~  0.  Eqs.  (22)  and  (29)  time 
scales  differ  by  a  (a^+n)  factor  resulting  in  a  1/100  ratio  in  Fig.  5 
(aj+n  =  0.1 ).  Steady  states  are  indicated  and  given  by  Eqs.  (24)  and 
(31 )  whereas  initial  linear  regimes  are  given  by  Eqs.  (23)  and  (30). 


4.  Results  and  discussion 


4.1.  Model  validation 


Two  problems  complicate  model  validation.  In  the  first  instance, 
ohmic  loss  measurements  do  not  represent  the  actual  ionomer  con¬ 
ductivity  [11].  During  steady  state  operation,  the  inactive  foreign 
cation  flux  is  equal  to  zero.  However,  ohmic  measurements  per¬ 
formed  with  alternating  currents  or  potentials  affect  all  cations 
present  in  the  ionomer.  Therefore,  the  resulting  measurement  val¬ 
ues  include  a  contribution  from  foreign  cations.  If  direct  current 
or  potential  methods  are  used  instead  to  ensure  establishment 
of  a  steady  state  foreign  cation  concentration  profile  (zero  flux), 
proton  depletion  occurs  at  one  of  the  electrodes  leading  to  ther¬ 
modynamic,  kinetic  and  mass  transport  losses  thus  complicating 
separation  of  the  ionomer  ohmic  loss  contribution. 

The  absence  of  complete  and  extensive  cell  performance  data 
sets  constitutes  the  second  validation  problem.  Only  few  data  sets 
are  available  in  the  form  of  current,  potential  or  resistance  time 
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Fig.  6.  Effect  of  NH3  contamination  on  transient  membrane  conductivity.  Model 
curves  result  from  fitting  to  Eq.  (18).  Parameter  values  are  listed  in  Table  2. 


Fig.  7.  Effect  of  gas  phase  relative  humidity  on  transient  membrane  conductivity. 
Model  curves  result  from  fitting  to  Eq.  (18).  Parameter  values  are  listed  in  Table  2. 


series  for  a  restricted  range  of  foreign  cations  (NH3/NH4+,  Na+, 
Co2+,  Fe3+,  Al3+,  Pt  (unspecified  valence)  [51-57]).  More  specifically, 
data  sets  are  generally  not  obtained  under  potentiostatic  control 
or  until  steady  states  are  achieved  (including  recovery).  Therefore, 
the  model  limit  of  applicability  cannot  be  quantitatively  defined 
because  relevant  experimental  data  sets  are  not  available  and,  a 
comparison  between  model  predictions  and  experimental  data, 
cannot  be  completed.  It  is  noted  that  cell  performance  losses  are 
significant  as  demonstrated  by  NaCl  poisoning  leading  to  a  current 
density  loss  of  33%  at  0.6  V  [55].  For  ammonia  (NH4+),  cell  per¬ 
formance  losses  were  concentration  dependent  and  amounted  to 
current  density  decreases  at  0.5  V  of  25  and  50%  for  respectively  1 3 
and  130  ppm  [55]. 

Recently,  poisoning  experiments  were  performed  with  NFI3,  an 
important  hydrogen  fuel  contaminant  [58],  and  a  membrane  rather 
than  a  complete  fuel  cell  membrane/electrode  assembly  [30,50,59]. 
Ohmic  measurements  were  completed  in  absence  of  reactants  and 
with  alternating  currents  or  potentials  (uniform  concentration  pro¬ 
files).  Results  are  thus  reflective  of  the  ionomer  state  (composition 
is  the  focus  rather  than  a  measurement  of  the  ohmic  loss  during  fuel 
cell  operation).  These  data  sets  offer  an  opportunity  to  partly  test 
the  present  ohmic  loss  contamination  model  without  interference 
from  other  loss  mechanisms.  Flowever,  recovery  cannot  be  investi¬ 
gated  in  this  manner  unless  liquid  water  is  circulated  and  contacts 
the  membrane. 

Figs.  6-8  illustrate  the  effects  of  NH3  concentration,  relative 
humidity  (RH)  and  temperature  on  membrane  conductivity  as  a 
function  of  time.  Eq.  (18)  with  4  parameters  (a,  b,  a^nt  qdjdf/qDj) 
was  used  and  fitted  to  experimental  data.  Parameters  a  and  b 
are  already  known  and  are  respectively  derived  from  the  initial 
(yx+n  =  0,  <7  =  a  +  b)  and  steady  state  (yx+n  =  1,  <7  =  a)  conductivity 
values.  The  two  other  parameters  were  obtained  by  least  square 
fitting  (SigmaPlot  version  10,  Marquardt-Levenberg  algorithm). 


Fig.  8.  Effect  of  temperature  on  transient  membrane  conductivity.  Model  curves 
result  from  fitting  to  Eq.  (18).  Parameter  values  are  listed  in  Table  2. 


Fig.  9.  Model  parameter  cgDgICididg  (Table  2)  as  a  function  of  cg.  A  significant  linear 
correlation  is  also  displayed. 


Table  2  lists  all  parameter  values  determined  using  Figs.  6-8  data. 
It  is  noted  that  correlation  coefficients  generally  exceed  0.97  indi¬ 
cating  that  the  model  accurately  reproduces  all  data  sets.  Clearly, 
NH3  diffusion  in  the  gas  phase  does  not  appear  to  correspond  to 
the  model’s  physical  picture  of  a  foreign  cation  moving  into  water 
(Fig.  la).  This  apparent  discrepancy  is  rationalized  by  pointing  out 
that  the  NFI3  species  still  moves  predominantly  by  diffusion  in  the 
gas  phase  boundary  layer  (absence  of  convection  and  migration). 
Thus,  in  Eqs.  (7)  and  (8),  liquid  phase  values  can  be  replaced  by 
their  gas  phase  equivalents.  Flere,  q  »  cg  and  d;  <  dg,  but  D;  <  Dg. 
However,  contamination  by  NH3  is  gas  phase  controlled  (Eq.  (8) 
applies)  because  qD/dg/CgDgdj » 1  and  5  +  2a^+n  >  1.  At  the  gas 
phase/ionomer  interface,  the  following  fast  reaction  [60]  ensures 
that  a  NH4+  concentration  is  established  effectively  creating  an 
effect  similar  to  ion  exchange  (a  foreign  cation  penetrates  the 
ionomer  whereas  the  proton  is  removed  from  the  ionomer,  Fig.  1  a): 

NH3  +  H+^NH4+  (32) 

The  reverse  reaction  is  not  considered  because  recovery  was  not 
observed  in  the  absence  of  ammonia  even  after  a  period  of  one  day 
[30].  A  final  confirmation  is  provided  by  plotting  the  CgDg/qd/dg 
parameter  (Table  2)  as  a  function  of  cg  (Fig.  9).  As  expected,  a  sig¬ 
nificant  linear  correlation  is  obtained.  Eq.  (32)  also  takes  place  in 
liquid  water  which  means  that  in  an  application  (in  the  presence 
of  a  membrane/electrode  assembly  rather  than  a  membrane),  the 
ammonia  species  behavior  is  closer  to  a  foreign  cation  (ammo¬ 
nium). 

4.2.  Model  parameters  determination 

All  Eq.  (22)  parameters  (a,  b ,  a^+n,  qdjdf/qDj)  are  easily  deter¬ 
mined  if  the  initial  ionomer  proton  form  conductivity  is  known  (Eq. 
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Table  1 

Ion  exchange  isotherm  parameters3  summary  (Fig.  3). 


Foreign  cation 

cxH+ 

wx+n 

r2 

nh4+ 

0.427 

0.988 

Na+ 

1.25 

0.998 

Ca2+ 

0.0933 

0.981 

Fe3+ 

0.00106 

0.965 

3  Nation  417  data  [31]. 


(16)): 

Est^=a  +  b  (33) 


The  limiting  dimensionless  current  density  during  contamina¬ 
tion  (Eq.  (24))  provides  the  experimental  value  required  to  isolate 
a  and  b  parameters: 


Est2  = 


a 

a  +  b 


(34) 


Eqs.  (33)  and  (34)  leads  to: 


a  =  Est^Est2 

(35) 

b  =  £sti(l  -  Est2) 

(36) 

Eqs.  (23)  and  (30)  have  a  similar  slope  providing  independent 
estimates  for  C/d/df/c/D/  (estimates  can  be  used  to  verify  model  con¬ 
sistency).  Eqs.  (33)  and  (36)  lead  to: 


Est3  = 


b  cjDi 

a  +  b  qdidi 


=  (1  -Est2) 


cfii 

C[d[d\ 


(37) 


Rearrangement  of  Eq.  (37)  leads  to: 


qdjd/  _  1  -  Est2 
C/D/  Est3 


(38) 


The  time  required  for  the  dimensionless  current  to  decrease  to 
a  value  equal  to  half  the  difference  between  initial  and  steady  state 
values  t0.5  (contamination)  represents  the  last  estimate: 


Est4  =  tp. 5 


(39) 


The  corresponding  dimensionless  current  is  using  Eq.  (24): 
i(to.s)  i(0)-i(oo)  1/  a  \ 

i(0)  2i(0)  2  V  a  +  b)  1  J 

Replacement  of  Eq.  (40)  in  Eq.  (22)  and  use  of  Eqs.  (35),  (36)  and 
(38)  leads  to: 


1  Est2-Est3Est4 


aH+  - 


9  + 


1  -Est2 


Ki-t! 


Est2 


Est2  y 


+  ^  + 


Est2 
1  -Est2 


(41) 


The  equivalent  derivation  for  the  recovery  case  using  Eqs.  (29) 
leads  to: 


VH+ 

xX+n 


ln(2  +  ll 


Est2 


Est2 


I  1  _  Est2 
^  2  1  -Est2 


Est2+Est3Est4 
1  ~Est2 


(42) 


Eq.  (42)  provides  a  second  estimate  that  can  be  used  to  verify 
model  consistency. 


4.3.  Contamination  mechanism 


It  was  indicated  in  Section  4.1  that  validation  data  sets  are 
incomplete.  However,  the  information  available  supports  the 
present  model.  For  instance,  partial  or  complete  performance 
recovery  was  noted  [51,53].  More  specifically,  experimental  data 
for  Fe3+,  representing  a  worse  case  scenario  with  «:  1  (higher 
charge  cations  are  strongly  preferred  by  the  ionomer  as  shown  in 
Fig.  3  and  Table  1  [26,31]),  are  reexamined.  A  membrane  in  the 
Fe3+  form  was  used  to  make  a  membrane/electrode  assembly  that 


was  conditioned  with  wet  reactant  streams  [14].  The  cell  voltage 
increased  from  0  V  until  cell  operation  was  possible  demonstrat¬ 
ing  ion  exchange  by  protons  from  water  dissociation  (Eq.  (1)).  A 
more  complete  model  validation  requires  additional  data  obtained 
under  specific  operating  conditions  (cell  voltage  control,  sequential 
contamination  and  recovery  transients  obtained  until  steady  states 
are  reached,  data  obtained  for  relevant  ranges  of  foreign  cation 
concentration  and  cell  potential,  etc.). 

A  model  library  is  systematically  and  incrementally  developed. 
As  the  current  density  is  increased,  leading  performance  losses 
change  from  kinetic  to  ohmic  and  mass  transport  losses  (Fig.  4).  The 
development  of  this  library  is  following  a  similar  path.  Initially,  con¬ 
taminants  leading  to  kinetic  losses  were  examined  [4,6]  followed 
by  ohmic  losses  (this  work).  Contaminants  creating  mass  transport 
losses  will  be  examined  after  the  library  is  populated  with  known  or 
envisaged  contaminant  mechanisms  leading  to  ohmic  losses.  Some 
of  the  existing  models  will  likely  need  a  revision  including  ionomer 
contamination  by  a  foreign  cation  to  take  account  of  thermody¬ 
namic,  kinetic  and  mass  transport  losses  associated  with  cations 
redistribution  within  the  ionomer  (this  work).  This  approach  was 
selected  to  manage  the  significant  amount  of  work  necessary  to 
cover  existing  ([4,6],  this  work)  and  envisaged  cases  by  separating 
it  into  more  easily  tackled  parts.  This  approach  is  also  expected  to 
facilitate  the  study  of  contaminant  mixtures  and  the  simultaneous 
presence  of  several  types  of  performance  losses. 

Table  3  summarizes  the  features  of  the  present  model  as  well 
as  other  model  versions  derived  for  other  cases  using  a  simi¬ 
lar  approach.  Most  of  these  cases  are  representative  of  known 
contamination  kinetics.  The  first  4  cases  are  related  to  catalyst  con¬ 
tamination  and  within  this  category,  3  cases  pertain  to  electroactive 
contaminants:  the  contaminant  X  adsorbs  and  reacts  leading  to  a 
product  P2  (such  as  CO  on  Pt  above  0.5  V  vs.  RHE  or,  NO  on  Pt  below 
0.6  V  vs.  RHE  or  above  0.9  V  vs.  RHE)  with  a  surface  reaction  or 
product  desorption  rate  determining  step  (rds),  the  contaminant  X 
adsorbs  but  only  reacts  at  segregated  sites  leading  to  an  irreversibly 
adsorbed  product  P2  (such  as  S02  or  H2S  on  Pt  between  0.05  and  1 V 
vs.  RHE).  The  4th  case  related  to  catalyst  contamination  represents 
an  electroinactive  contaminant  (such  as  CO  on  Pt  at  <0.5  V  vs.  RHE  or 
NO  on  Pt  between  0.6  and  0.9  V  vs.  RHE).  The  5th  case  pertains  to  the 
present  model  and  expands  for  the  first  time  the  model  library  into 
ionomer  contamination  kinetics  (a  significant  and  important  step). 
The  current  density  i  qualitatively  shows  a  similar  transient  behav¬ 
ior  for  all  cases  due  to  the  overwhelming  effect  of  the  much  more 
concentrated  reactant  species  with  the  exception  of  an  electroac¬ 
tive  contaminant  that  leads  to  an  irreversibly  adsorbed  product  P2. 
Thus,  a  distinction  between  all  cases  is  not  possible.  By  contrast, 
experiments  performed  without  a  reactant  lead  to  the  contami¬ 
nant  related  current  density  ix  or  ix+n  and  different  transients  for 
each  case  with  the  exception  of  the  present  membrane  contamina¬ 
tion  model.  Contamination  mechanism  identification  with  minimal 
testing  requirements  is  thus  significantly  facilitated.  Separation 
between  membrane  and  catalyst  (X  electroinactive)  contamination 
cases  is  conceivably  possible  using  additional  cell  resistance  mea¬ 
surements  [51  ].  An  accurate  measurement  is  not  necessary  (Section 
4.1)  but  rather  only  an  indication  that  the  membrane  resistance 
is  affected  to  discriminate  between  cases.  As  already  indicated  in 
Section  2,  cation  distribution  gradients  within  the  ionomer  prevent 
accurate  membrane  resistance  measurements  [7-13,61,62]. 

4.4.  Model  predictions 

It  was  reported  that  small  amounts  of  foreign  cations,  such  as 
Na+,  Cu2+,  Ni2+,  Fe3+,  significantly  impact  the  membrane  conduc¬ 
tivity  at  the  ~1 0-310  ppm  level  [63-65].  Therefore,  there  is  a  need 
to  prevent  the  gradual  replacement  of  ionomer  protons  by  foreign 
cations  because  procedures  to  maintain  such  a  low  contamina- 
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Fig.  10.  Effect  of  the  separation  factor  on  dimensionless  contamination  transients. 
a  =  0.01  £2_1  cm-1  and  b  =  0.09  cm-1 . 


cfd,c/(a"l 


Fig.  11.  Effect  of  the  ionomer  conductivity  on  dimensionless  contamination  tran¬ 
sients.  a  and  b  units  are  cm-1,  =0.1. 


tion  level  are  still  under  development.  The  ionomer  contamination 
model  suggests  operating  conditions  and  design  based  mitigation 
strategies.  Eq.  (22)  indicates  that  a  decrease  in  temperature  (D/) 
and  a  relative  decrease  in  foreign  cation  concentration  (q)  slow 
the  contamination  process.  Conversely,  an  increase  in  temperature 
and  a  relative  increase  in  proton  concentration  (relative  decrease  in 
foreign  cation  concentration)  accelerate  the  recovery  process  (Eq. 
(29)).  From  a  design  standpoint,  strategies  are  more  limited  because 
practical  ionomer  choices  (q)  and  their  thickness  (d2-)  are  restricted 
by  other  considerations.  For  instance,  lower  equivalent  weights 
and  thicker  materials  are  respectively  more  prone  to  water  uptake 
(induced  mechanical  stresses)  and  larger  ohmic  losses.  Similar  con¬ 
siderations  apply  to  the  gas  diffusion  layer  (d/)  with  a  negative 
impact  on  volumetric  power  density  and  material  cost. 

Eqs.  (22)  and  (29)  indicate  that  a,  b  and  parameters  influ¬ 
ence  transients  in  several  ways  because  they  appear  multiple  times. 
Thus  the  overall  effect  is  not  easily  discernible.  As  a  result,  the 
effect  of  all  these  parameters  on  both  contamination  and  recov¬ 
ery  transients  were  computed  for  a  relevant  range  of  parameter 
values.  Although  relative  humidity  is  not  specifically  included  in 
the  model  derivation,  model  predictions  based  on  ionomer  con¬ 
ductivities  resulting  from  exposure  to  different  relative  humidity 
environments  are  expected  to  hold  as  long  as  liquid  water  is  present 
in  the  flow  field  channel/gas  diffusion  electrode  and  contacts  the 
ionomer  because  the  physical  description  of  the  process  still  holds 
(Fig.  la).  Fig.  10  illustrates  the  effect  of  which  shows  that  a 
smaller  value  delays  contamination.  This  is  at  first  sight  counterin¬ 
tuitive  because  lower  separation  factor  values  are  associated  with 
higher  valence  foreign  cations  which  are  strongly  preferred  by  the 
ionomer  (Fig.  3).  However,  the  time  scale  q dj-d/a^/c/D/  (Eq.  (22))  is 
in  effect  largely  reduced  as  expected  thus  creating  a  larger  dimen¬ 
sionless  time.  It  is  also  noted  that  the  approach  to  a  steady  state 
value  is  sharper  as  the  separation  factor  value  decreases.  Fig.  1 1 
displays  the  effect  of  both  relative  humidity  and  foreign  cation. 
The  top  2  curves  correspond  to  an  ionomer  exposed  to  a  saturated 
stream  ( a  +  b  =  0.1  cm-1  [31  ])  whereas  the  bottom  2  curves  cor¬ 
respond  to  an  ionomer  exposed  to  a  50%  relative  humidity  stream 
(a  +  b  =  0.04  cm-1  [31]).  For  the  same  relative  humidity,  each 
curve  pair  shows  the  range  of  expected  ionomer  conductivity  in 
the  foreign  cation  form  (a  =  0.01  for  Na+,  Ca2+,  Fe3+  to  0.06  for 
NH4+  and  0.0004  for  Na+,  Ca2+,  Fe3+  to  0.003  for  NH4+  cm-1 
[31  ]).  It  is  observed  that  both  relative  humidity  and  foreign  cation 
do  not  influence  the  contamination  process  as  expected  because 
the  time  scale  qd/d/a^+n/c/D/  is  not  dependent  on  these  variables. 
It  is  also  observed  that  relative  humidity  has  a  large  impact  on 
the  dimensionless  current  density  steady  state  value.  Therefore, 
ionomer  contamination  by  foreign  cations  is  relatively  more  severe 


Fig.  12.  Effect  of  the  separation  factor  on  dimensionless  recovery  transients. 
a  =  0.01  £2_1  cm-1  and  b  =  0.09  cm-1 . 

under  low  relative  humidity  conditions  and  does  not  support  the 
development  of  ionomers  able  to  operate  under  such  conditions 
[66].  Fig.  1 2  shows  the  effect  of  ot^+n  on  recovery  transients  demon¬ 
strating  that  larger  values  delay  recovery  to  a  proton  form.  This 
counter-intuitive  result,  a  small  separation  factor  is  expected  to 
facilitate  ion  exchange  by  protons,  is  explained  by  considering  Eq. 
(29)  time  scale  cldId//c/D/G'11+n.  A  larger  separation  factor  leads  to 
a  smaller  time  scale  and  thus  to  a  larger  dimensionless  time  value. 
The  approach  to  a  steady  state  value  remains  smooth  for  all  sepa¬ 
ration  factor  values  contrasting  with  Fig.  10  results.  Fig.  13  depicts 
the  effect  of  both  relative  humidity  and  foreign  cation  in  a  similar 


cp,ah;j 

cp'd i 


Fig.  13.  Effect  of  the  ionomer  conductivity  on  dimensionless  recovery  transients,  a 
and  b  units  are  £2_1  cm-1,  a^_n  =0.1. 
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Table  2 

Contamination  model  parameters  summary. 


Validation  data 

a  cm"1) 

cm-1) 

OfH+ 

“x+n 

CididilqDi/h 

r2 

Fig.  6 

5  ppm  NH3 

0.00261 

0.0270 

0.0617 

27.2 

0.990 

10  ppm  NH3 

0.00234 

0.0253 

0.0394 

20.7 

0.994 

20  ppm  NH3 

0.00255 

0.0265 

0.131 

5.75 

0.997 

30  ppm  NH3 

0.00225 

0.0244 

0.281 

2.30 

0.983 

Fig.  7 

30%  RH 

0.00023 

0.0113 

0.752 

1.44 

0.973 

50%  RH 

0.00248 

0.0269 

0.217 

7.08 

0.856 

65%  RH 

0.00737 

0.0465 

0.210 

8.87 

0.974 

80%  RH 

0.0205 

0.0763 

0.0761 

10.3 

0.916 

90%  RH 

0.0337 

0.0739 

0.0805 

9.46 

0.969 

100%  RH 

0.0552 

0.131 

0.0286 

16.3 

0.990 

Fig.  8 

60  °C 

0.00193 

0.0247 

0.211 

4.91 

0.971 

80°C 

0.00254 

0.0265 

0.108 

6.03 

0.993 

90  °C 

0.00358 

0.0297 

0.157 

3.72 

0.973 

Table  3 

Summary  of  model  behavior  for  different  fuel  cell  contamination  kinetics. 

Measured  variable  for  step  changes  in  contaminant  concentration 

i  (cr  +  0) 

ix  or  ix+n  (cR  =  0) 

Catalyst 

contamination 

kinetics 


X  electroactive 
leading  to  a 
product  P2 


X  electroinactive 
[4] 


X  reaction  rds  [4] 


P2  desorption  rds  [4] 


Irreversible  P2 
adsorption  [6] 


©  © 

It,  t 

©-© 

Catalyst  surface 

©  © 

J  rds\ 

©~»© 

Catalyst  surface 

©  © 

it  l 

©  ©-© 

Segregated  catalyst  surface 

© 

It 

© 

Catalyst  surface 


t 


t 


'x 

0 - 


Membrane  X+n  electroinactive 

contamination  (this  work) 

kinetics 


Liquid 

water 


lonomer 


X+n 


n  H+ 


0 


t 


fashion  as  Fig.  11.  It  is  observed  that  in  each  case  full  recovery  is 
achieved  with  a  time  scale  independent  of  relative  humidity  and 
foreign  cation  despite  a  large  variation  in  initial  dimensionless  cur¬ 
rent  density. 

Recently,  low  levels  of  foreign  divalent  or  trivalent  cations  were 
proposed  to  mitigate  ionomer  degradation  (radical  scavengers) 
[67-69].  Experimental  data  show  that  these  cations  are  retained 
by  the  ionomer  for  a  period  of  time  exceeding  6000  h  [67,69].  The 
present  model  properly  predicts  the  long  residence  time  of  a  higher 
valence  cation  in  an  ionomer.  For  example,  Fig.  12  shows  a  dimen¬ 
sionless  time  of  about  4  for  recovery  in  the  presence  of  a  trivalent 
cation  (a£+n  ~0.001  in  Table  1 ).  The  combination  of  this  value  with  a 


~0.001  separation  factor  for  Fe3+  (Table  1 )  and  a  value  of  C/d/dj/c/D/ 
of  ~l-30h  (Table  2)  leads  to  a  recovery  time  varying  from  4000  to 
1 20,000  h  which  is  in  good  agreement  with  the  observed  6000+ h 
time  scale. 

5.  Conclusion 

A  PEMFC  contamination  model  was  derived  for  the  case  of 
a  foreign  cation  displacing  the  ionomer  proton.  The  model  was 
derived  by  taking  account  of  only  key  processes  to  obtain  analyt¬ 
ical  solutions  and  facilitate  a  widespread  implementation.  Model 
validation  indicates  that  the  simplified  approach  derivation  is 
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accurate  and  represents  a  benchmark  for  a  more  complicated 
model  that  includes  migration  effects.  Additional  fuel  cell  val¬ 
idation  data  obtained  at  low  and  high  current  densities  (cell 
voltage  control,  sequential  contamination  and  recovery  transients 
obtained  until  steady  states  are  reached)  are  desirable  including 
the  effect  of  foreign  cation  composition  and  concentration,  etc. 
This  is  especially  important  considering  many  contaminants  and 
their  effects  have  not  yet  been  identified  or  defined  including 
cations  that  originate  from  system  materials  or  are  produced  from 
reactions  involving  airborne  species  [2].  These  cation  sources  are 
currently  under  investigation  within  the  United  States  Department 
of  Energy  program.  Other  model  extensions  need  to  be  considered 
including  foreign  cations  that  are  reducible  and  that  impact  oxygen 
reduction  kinetics  (Cu2+  and  Co2+  [12],  Ru+n  [70]).  The  effect  of  the 
foreign  cation  uneven  distribution  within  the  ionomer  affecting 
thermodynamic,  kinetic  and  mass  transport  losses  also  need  to 
be  taken  into  account  [7-13].  From  that  standpoint,  the  negative 
impact  of  foreign  cations  in  the  ionomer  on  reactant  permeability 
[71-75]  has  not  yet  been  considered.  Model  predictions  suggest 
the  use  of  mitigation  strategies  and  a  reconsideration  of  low  water 
requirement  ionomer  development  programs. 
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